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1 Introduction 

The physics of the 20th century is the physics of local gauge invariance. Lattice 
gauge theories present an extreme limit, being systems with exact internal 
symmetries but no (or very restricted) space time symmetries. These lectures 
are intended to be an introduction to the lattice approach to quantum field 
theories in general and to QCD in particulartl. The use of a lattice regulator is 
a standard technique for studying strong coupling when space-time symmetries 
can be regarded as secondary. Many techniques of lattice have found a home 
or at least a word in the vocabulary of string. Some of these techniques might 
be useful to you. Lattice methods have been used to study random surfaces in 
the context of models for quantum gravity, but I will not discuss that subject 
here. 

Lattice gauge theory was invented by Wilson □ in 1974. For its first five 
years or so it was a subjeci studied via analytic calculations. Then in 1979 
Creutz, Jacobs, and RebbiH performed the first numerical studies of a lattice 
gauge theory using Monte Carlo simulation, and over the next few years the 
field was transformed into one in which numerical methods combined with 
analytic methods to study nonperturbative behavior. By the late 1980's the 
field was dominated by very professional large scale simulations of QCD on 
supercomputers. In the last few years some members of the lattice community 
have been developing techniques for reducing the size of numerical simulations. 
These techniques rely heavily on modern work on the renormalization group 
and effective actions. Some of them actually work: some lattice calculations 
can now be done on work stations. Some of the techniques may have a wider 
range of applicability than lattice QCD. 

Most of the theoretical ideas that this audience might find interesting are 
from the pre-computer days of lattice gauge theory. I will begin by describing 
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how gauge symmetry is realized on the lattice and then show how one studies 
confinement in the strong coupling limit. Next, I will describe strong- weak 
coupling duality in the context of lattice models. (This is the only place where I 
know how to explicitly perform a duality transformation on the field variables.) 

1 will next describe how lattice methods are used to obtain qualitative insight 
into the behavior of various model systems, without using the computer. Then 
the subject will change to QCD. I will describe how one performs Monte Carlo 
simulations, and show the state of the art in lattice calculations of hadron 
spectroscopy. In one lecture I will talk about the connection between bare 
actions, effective actions, and scaling, and how this connection is being used 
to construct "improved actions" for QCD. I will present examples to show the 
extent to which they actually do improve. Two case studies: glueballs, and 
the calculation of as{Mz) from the lattice, will conclude the review. 

2 Gauge Field Basics 

2. 1 Beginnings 

The lattice is a cutoff which regularizes the ultraviolet divergences of quantum 
field theories. As with any regulator, it must be removed after renormalization. 
Contact with experiment only exists in the continuum limit, when the lattice 
spacing is taken to zero. 

We are drawn to lattice methods by our desire to study nonperturba- 
tive phenomena. Older regularization schemes are tied closely to perturbative 
expansions: one calculates a process to some order in a coupling constant; 
divergences are removed order by order in perturbation theory. The lattice, 
however, is a nonperturbative cutoff. Before a calculation begins, all wave- 
lengths less than a lattice spacing are removed. Generally one cannot carry 
out analytical studies of a field theory for physically interesting parameter val- 
ues. However, lattice techniques lend themselves naturally to implementation 
on digital computers, and one can perform more-or-less realistic simulations of 
quantum field theories, revealing their nonperturbative struc;ture, on present 
day computers. I think it is fair to say that little of the quantitative results 
about QCD which have been obtained in the last decade, could have been 
gotten without the use of numerical methods. 

On the lattice we sacrifice Lorentz invariance but preserve all internal 
symmetries, including local gauge invariance. This preservation is important 
for nonperturbative physics. For example, gauge invariance is a property of 
the continuum theory which is nonperturbative, so maintaining it as we pass to 
the lattice means that all of its consequences (including current conservation 



2 



and renormalizability) will be preserved. 

It is very easy to write down an action for scalar fields regulated by a 
lattice. One just replaces the space-time coordinate a;^ by a set of integers 
(xp = an^, where a is the lattice spacing). Field variables (j){x) are defined on 
sites 4>{xn) = 4>rn The action, an integral over the Lagrangian, is replaced by 
a sum over sites 



y"d^x/:^a4^£(0„)- (1) 



and the generating functional for Euclidean Green's functions is replaced by 
an ordinary integral over the lattice fields 

{\{d4>n)e-^' ■ (2) 

n 

Gauge fields are a little more complicated. They carry a space-time index /i in 
addition to an internal symmetry index a and are associated with a 

path in space Xfi{s): a particle traversing a contour in space picks up a phase 
factor 



P(exp ig j^dXf^Af_,)ilj (3) 
= U{s)^{x). (4) 



P is a path-ordering factor analogous to the time-ordering operator in ordinary 
quantum mechanics. Under a gauge transformation g, U{s) is rotated at each 
end: 

U{.s) ^ g-\x^is))U{s)g{x^m. (5) 

These considerations led Wilson i to formulate gauge fields on a space-time 
lattice, as follows: 

The fundamental variables are elements of the gauge group G which live 
on the links of a D-dimensional lattice, connecting x and x + fi: U^{x), with 

C/^(n)=exp(*5aT''A;:(n)) (6) 

for SU{N). {g is the coupling, a the lattice spacing, the vector potential, 
and T" is a group generator). 

Under a gauge transformation link variables transform as 

U^{x)^V{x)U^{x)V{x + t,y (7) 

and site variables as 

^(x) -> V{x)il;{x) (8) 
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so the only gauge invariant operators we can use as order parameters are matter 
fields connected by oriented "strings" of U's (Fig. la) 

iPixi)Uf,{x)U^{x + fi)...^{x2) (9) 

or closed oriented loops of U's (Fig. lb) 

Tr . . . U^ix)U^ix + fi) . . . ^ Tr . . . U^{x)V^ ix + fi)Vix + fi)U^{x + fi) .... (10) 

An action is specified by recalling that the classical Yang-Mills action in- 
volves the curl of , F^^i, . Thus a lattice action ought to involve a product of 
C/^'s around some closed contour. In fact, there is enormous arbitrariness at 
this point. We are trying to write down a bare action. So far, the only require- 
ment we want to impose is gauge invariance, and that will be automatically 
satisfied for actions built of powers of traces of U's around closed loops, with 
arbitrary coupling constants. If we assume that the gauge fields are smooth, 
we can expand the link variables in a power series in gaA'^s. For almost any 
closed loop, the leading term in the expansion will be proportional to F^^. We 
might want our action to have the same normalization as the continuum action. 
This would provide one constraint among the lattice coupling constants. 

The simplest contour has a perimeter of four links. In SU {N) 

IN 

/35=— ^^ReTr {l - U^{n)U,{n + fi)Ul{n + v)Ul{n)). (11) 

This action is called the "plaquette action" or the "Wilson action" after its 
inventor. _ 

Note Elitzur's theorem: □ only gauge invariant quantities have nonzero 
expectation values. 

Let us see how this action reduces to the standard continuum action. Spe- 
cializing to the U(l) gauge group, 

5 = 4 51 51 {'^-eM^9a[Mn)+A^{n+fi)~Af,{x + i))-A^{n)])). (12) 

The naive continuum limit is taken by assuming that the lattice spacing a is 
small, and Taylor expanding 

Af,{n + u) = Af,{n) + ad^Af,{n) + . . . (13) 

so the action becomes 

S = ^ E„ E^>. 1 - Re {eM^ga[a{^.A,, - d^^A,) -f 0{o?)\)) (14) 
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(a) 



(b) 



Figure 1: Gauge invariant obscrvablcs arc cither (a) ordered chains ("strings") of links 
connecting quarks and antiquarks or (b) closed loops of link variables. 
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(15) 
(16) 



transforming the sum on sites back to an integral. 

The best way people have found to check for quark confinement is to ask 
the following question: In a world in which there are no light quarks, what is 
the potential V{R) between a heavy qq pair? If the limit as R goes to infinity 
of V{R) is infinite, we have confinement, if not, quarks are not confined. V{R) 
can be computed by considering the partition function in Euclidean space for 
gauge fields in the presence of an external current distribution: 



(17) 



(for a non-Abelian gauge group insert color sums J°A°.) If represents a 
point particle moving along a world line, it is a ^-function on that world line 
(parameterized by £^): 



i j Jij,Aij,d^x = i j> A^dlfj,. 



(18) 
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R 1 2 

(a) (b) 



Figure 2: An R X T Wilson loop (a) and the pattern of tiling which occurs when it is 

evaluated in strong coupling (b). 

Gauge invariance, which imphes current conservation, says that the current 
Une cannot end, and so the loop is a closed loop. Let's make it rectangular 
for simplicity, extending a distance R in some spatial direction and a distance 
T in the Euclidean time direction. This is the famous Wilson loop. (See Fig. 
2a.) 

Zj describes the loop immersed in a sea of gluons. We can think of it as 
the partition function for a qq pair which is created at i = 0, pulled apart a 
distance R, and then allowed to annihilate at a time T later. 

To find the pair's energy, differentiate ln(Z) with respect to T (physicallv, 
we are measuring the response of the system as we stretch the loop a bit). U 
But this is not quite right. This calculation will include the vacuum energy 
of the gluons, which diverges. What we really want is the change in energy of 
the system when we include the pair and when it is not present: 

dT '■^ J[dU]e-I^S 

= -^ln(M^)j=o 
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(19) 
(20) 

(21) 
(22) 



where (W) is the expectation value of the Wilson loop in the background gluon 
field. Thus the behavior of (W) will tell us whether we have confinement or 
not. 

For example, if (W) ~ exp(— m(2i? + 2T)) (a so-called "perimeter law") 
then E = 2m, and m can be interpreted as the quark mass. The energy 
of the pair is independent of R so that the quarks can separate arbitrarily 
far apart. However, if (W) ~ exp(— cri?T) (an "area law"), E = aR, and 
quarks are confined with a linear confinement potential. The parameter a 
is called the "string tension" . In general. In W is complicated: —ln{W) = 
(jRT + 2m{R + T) + constant + T/i? + . . . Nevertheless, any area law term 
will dominate at large R and T and the theory will confine. Notice that this 
confinement test will fail in the presence of light dynamical fermions. As 
we separate the heavy QQ pair, at some point it will become energetically 
favorable to pop a light qq pair out of the vacuum, so that we are separating 
two color singlet mesons. The Wilson loop will show a perimeter law behavior. 

Now; let us give an explicit demonstration of confinement for U(l) gauge 
theory, cl The Haar measure for U{\) is {6^ — go-A^) 

Z = j [dU] exp{~pS) (24) 

where S is given by Eq. |ll|. 

Strong coupling corresponds to small /?. Now we write down two trivial 
(but useful) results: 

" ^^^i9+4>) ^ (25) 

i(e+0l)g-i(9+02) ^ gi(0l-02) (26) 



and 



/' 

J — TT 



de 

2^' 



For small p we may expand the exponential in Z as a power series in f3: 

/Hf) 
n (l + /3(Tr[/p + h.c. + ...) (27) 



links plaquettes 

l + OiP^) (28) 



since all linear terms vanish. 
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Let us compute the expectation value of an L x T Wilson loop (Fig. 2a): 
There are LT plaquettes enclosed by it. 

{wmn^^JH^cM^T.^) n a+pcos{j2e) + ...) (29) 

WL plaquettes 

It's easy to convince one's self that W is proportional to P'"'^ . Pick one link 
on the boundary to integrate. We need a term in exp(— 5') to contribute and 
cancel the phase of the boundary link appearing in the Wilson loop (link 1 of 
Fig. 2b). Now we simultaneously need another 0(/3) term to cancel the phases 
of link 2 (of Fig. 2b). Only if we "tile" the whole loop as shown do we get a 
nonzero expectation value; each plaquette contributes a factor /3 so 

{W{L,T)) (30) 

Eqq = UL (31) 

cr =-ln^ (32) 

This is confinement by linear potential. The calculation can be generalized 
to all gauge groups. In the strong coupling limit, any lattice-gauge theory 
shows confinement. This calculation was first done by Wilson tl and was the 
first demonstration of quark confinement in a field theory in more than two 
space-time dimensions, albeit in an unphysical limit. This result had a lot to 
do with the early interest in lattice gauge theory: one has an approximation 
in which confinement is automatic. 

2.2 Taking The Continuum Limit, and Producing a Number in MeV 

When we define a theory on a lattice the lattice spacing a is an ultraviolet cutoff 
and all the coupling constants in the action are the bare couplings defined with 
respect it. When we take a to zero we must also specify how g(a) behaves. 
The proper continuum limit comes when we take a to zero holding physical 
quantities fixed, not when we take a to zero holding the couplings fixed. 

On the lattice, if all quark masses are set to zero, the only dimensionful pa- 
rameter is the lattice spacing, so all masses scale like 1/a. Said differently, one 
computes the dimensionless combination am{a). One can determine the lat- 
tice spacing by fixing one mass from experiment. Then all other dimensionful 
quantities can be predicted. 

Now imagine computing some masses at several values of the lattice spac- 
ing. (Pick several values of the bare parameters at random and calculate 
masses for each set of couplings.) Our calculated mass ratios will depend on 
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the lattice cutoff. The typical behavior will look like 



(ami(a))/(am2(a)) = mi(0)/TO2(0) + O(mia) + 0{{miaf) + ... (33) 

The leading term does not depend on the value of the UV cutoff, while the 
other terms do. The goal of a lattice calculation (like the goal of almost any 
calculation in quantum field theory) is to discover the value of some physical 
observable as the UV cutoff is taken to be very large, so the physics is in 
the first term. Everything else is an artifact of the calculation. We say that 
a calculation "scales" if the a— dependent terms in Eq. ^ are zero or small 
enough that one can extrapolate to a = 0, and gcnerically refer to all the 
a— dependent terms as "scale violations." 

We can imagine expressing each dimensionless combination am{a) as some 
function of the bare coupling(s) {g{a)}, am = f{{g(a))}. As a —> we must 
tune the set of coupHngs {g{a)} so 

lim —f({g(a)}) constant. (34) 
a— >o a 

From the point of view of the lattice theory, we must tune {g} so that correla- 
tion lengths 1/ma diverge. This will occur only at the locations of second (or 
higher) order phase transitions in the lattice theory. 
Recall that the /3-function is defined by 

dgja) dg{a) 

(There is actually one equation for each coupling constant in the set. A is 
a dimensional parameter introduced to make the argument of the logarithm 
dimensionless.) At a critical point /3((?c) = 0. Thus the continuum limit is the 
limit 

lim{5(a)} ^ (36) 

Continuum QCD is a theory with one dimensionless coupling constant. In 
QCD the fixed point is = so we must tune the coupling to vanish as a 
goes to zero. 

Pushing this a little further, the two-loop /3-function is prescription inde- 
pendent, 

(3{g)=-bi9'' + h2g\ (37) 

and so if we think that the lattice theory is reproducing the continuum, and 
if we think that the coupling constant is small enough that the two-loop beta- 
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function is correct, we might want to observe perturbative scaling, or "asymp- 
totic scaling", m/A fixed, or a varying with g as 



aA= ( 



1 



b2/(2fc?) 



exp( 



1 



)■ 



(38) 



9- 



big' (a) 



Please note that asymptotic scaling is not scaling. Scaling means that di- 
mensionless ratios of physical observables do not depend on the cutoff. Asymp- 
totic scaling involves perturbation theory and the definition of coupling con- 
stants. One can have one without the other. (In fact, one can always define a 
coupling constant so that one quantity shows asymptotic scaling.) 

And this is not all. There are actually two parts to the problem of produc- 
ing a number to compare with experiment. One must first see scaling. Then 
one needs to set the scale by taking some experimental number as input. A 
complication that you may not have thought of is that the theory we simulate 
on the computer is different from the real world. For example, a commonly 
used approximation is called the "quenched approximation" : one neglects vir- 
tual quark loops, but includes valence quarks in the calculation. The pion 
propagator is the propagator of a qq pair, appropriately coupled, moving in a 
background of gluons. This theory almost certainly does not have the same 
spectrum as QCD with six flavors of dynamical quarks with their appropriate 
masses. (In fact, an open question in the lattice community is, what is the 
accuracy of quenched approximation.) Using one mass to set the scale from 
one of these approximations to the real world might not give a prediction for 
another mass which agrees with experiment. We will see examples where this 
is important. 

3 Abelian Duality with a Lattice Regulator 

"Every theorist knows" that a phase with electrical confinement (area law for 
Wilson loops, etc.) is has a dual description in terms of magnetic variables 
which are screened, in analogy with the Meissner effect in a superconductor. 
There are many sources for this intuition. The one that people in my gener- 
ation point to are lattice calculations in Abelian gauge and spin models. In 
these models one can explicitly perform the duality transformation. In some 
models (notably the d = 2 U(l) spin model and the c? = 3 gauge model) one 
can solve the dynamics sufficiently well to understand the role of magnetic vari- 
ables: the Kosterlitz-Thouless vortex-unbinding transition in the former case, 
and PolviikQii's demonstration of confinement in the latter. Classic references 
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Let's begin with a. d = 2 U{1) (or 0(2)) spin model. On each site of the 

lattice there is a unit- length spin s(x), which can be labeled by an angle (the 
orientation of the spin) 4>{x). We assume nearest neighbor spins interact via 
H = J2xi^ ~ ^(*) ■ ^(^ + A)) = J2x,^ii^ ~ cos{(j){x) - (j){x + /*))). The partition 
function is 

Z= [d(t>{x)]exp{-pH). (39) 
Let's perform the Fourier transform 

oo 

exp ( - /3(1 - cos(0))) = J2 ^'"'fnW (40) 

— CSO 

with 

urn - ^e-""'^^''^ (41) 

to rewrite the partition function as 

Z = J[d(j)]'[[ e'X;"''(^)«^+/i)-0(a:))g-E,,^"M(^)V(2/3)^ |-42) 

We can now integrate over the spin variables, which constrains the nj[^s: 

j rf(/>(a;)e''^(^)["''(^-'')-"''(^)l = (5v^„^,o. (43) 

We can solve the constraint by introducing 

n^{x) = e^^V^n(a;) (44) 
which gives a partition function 

Z = n E exp(-E^), (45) 

X n{x) = -oo r,iJ, ^ 

but at large (3 the sums converge slowly. The Poisson resummation formula 



or 



J29il)=J2 fdM^y'"''"'^ (47) 
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allows us to rewrite the partition function as 



Z = l[lmx)] exp(-^(^|^)exp(27rz^m(r)0(r)) (48) 

or 

Z = ^m=0^m5^0 (49) 

where 

Zrn=o = I m^)] exp(- E ^^1^) (50) 

is a free massless field theory (these are "spin waves" in condensed-matter 
jargon), and the second term is a "vortex component," 

oo 

Z„^o = W ^ exp ( - 2Tr^ [3m{x)G{x - x')m{x')) (51) 

X m{x) — ~oo 

where 

V^G^Six^x'). (52) 

That is, we have rewritten the original partition function in terms of a set of 
variables which live on the sites, take on integer values (0, ±1, ±2, . . .), and 
interact by long-range Coulomb interactions. These variables are called (de- 
pending on your dialect) vortices, defects, solitons, monopoles, or instantons. 
Clearly, the /3 — > cx) phase has few vortices, and the low /3 phase will have 
many. This model has a^ phase transition at intermediate /3 mediated by the 
condensation of vorticesQ. 

The algebra is similar in other dimensions. For d ~ 3 the constraint is 

V^n^ = ^ = ep^AV^nA(x) (53) 

and if we impose a further "gauge constraint" 

VAriA = (54) 

we arrive at the partition function for a system of interacting, conserved, ori- 
ented "vortex strings" 

oo 

Z {string) =Y\_ ^ exp (-27r^/3 ^ ■rnf^{x)Gf^,i,{x-x')m^{x'))\^^rn^=o- 

JTif, (x) = — CXD x,x',fi,v 

(55) 
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Again, 



(56) 



and the oriented strings have a Coulomb interaction. 
U(l) gauge theories with the Wilson action 

cos(6l^^ (x) ) = 008(6*^ (x) +e^{x + ti)- 0^, {x + v)- 9^ (x) ) (57) 

can be treated similarly. The Fourier transform introduces a constraint 

exp(iZ^^6i^^) ^ (5v^i^„,o- (58) 

In three dimensions the constraint is solved by setting /^j^ = ^fj,v\d\l. We 
have a system of interacting (charged) point defects: monopoles or instantons. 
In four dimensions the constraint is solved by setting Z^i, = e^uxOxla with 
dcla = 0. This is a system of interacting monopole world lines. For more 
details about "what comes next," seel3ll3iy. 



3.1 Polyakov's Solution oj 3-d U(l) Gauge Theory 

Polyakov showed that the d = 3 U(l) gauge theory is permanently confining 
for all values of the gauge coupling. The calculation begins by considering the 
grand partition function for a gas of monopoles, where we restrict the possible 
charges to ±1. (In keeping with the spirit of the original demonstrations, 
all my numerical factors are probably wrong.) The monopole fugacity C, = 
exp(— 1/e^). The partition function is 

j-N f 

z-T. d^^---'^^^ ^M-P-" E (59) 



or 



Z= [[dx]eM~l I dx{Vx?)Y: E ^rf^i...d^^e'S^^^^^^^(^^\ (60) 

H J AT „ II 



N q=±l 



We can sum over the charges on a site J dx{e^^ + e '-^) and exponentiate 

to produce a sine-Gordon Lagrangian 

Z = J [dx]exp(-^ J dx[{Vxf -M^cosx]) (61) 
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where Iv'P — (3C,/ii. Now let's introduce a Wilson loop into the calculation. 
We reverse the derivation of the Wilson loop of the last lecture to write 

W{C)= eyipif^A-dl (62) 

= exp i J d'^xJext ■ A (63) 

= expiX;r<7(?')»7('') (64) 

where rj is the magnetic scalar potential arising from the external current loop, 
V^ry = V • Bext = P- We repeat the derivation including the external scalar 
potential and find 

Z{ti) - j [dx] exp(-^ j dx[{V{x - v)? " M'cosx])- (65) 

We solve the integral by steepest descent, and find 

Z(7?) = exp ( - ^ y dx[(V(xc/ - v)f - cos Xd]) (66) 

where 

^^Xci= sin Xci + ^^V (67) 
What is the solution? x only near the sheet and x — exp(— M|z|) moving 
transversely away from the sheet, so Z{ri) ~ exp(— 7 x area) where 7 ~ tt/ (3 x 
X 1/M ~ exp(— 1/e^). This is the signal of linear confinement. Hence 
the intuition: electrical confinement is due to the free monopoles which screen 
external magnetic fields. 



4 Simple Lattice Intuition 

Let's finally consider some variations on these models. Zn models are like 
f/(l) models, except the variables are restricted to 6 = 2Trj/N, j — 1 ... TV. 
They can have three phases: a low-beta phase with many defects, a high- 
beta magnetically-ordered or "Higgs" phase, and possibly a spin wave phase 
at intermediate coupling. The generic phase diagram is shown in Fig. ^ It is 
easy to estimate the location of the magnetic ordering transition: Fluctuations 
in go like (Sd)"^ ~ {9'^) ~ 1/(3, but if this 59 < \/N the spin cannot fluctuate 
from one allowed value to another one, and the system freezes. This happens 
at /3c ^ N^. 

While the spin-wave phase is "just free field theory," in two dimensions 
the free behavior of the angular variable is 

(<^,0,)^ilog(^). (68) 
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Figure 3: Phase diagram for two-dimensional Z[N) spin models. 



However, the spin-spin correlator is 

{Si-Si)= (e'W--<^.)) (69) 
= |r,,|(i//3) (70) 

whose exponent varies continuously with the coupling. The continuous vari- 
ation in a critical index with an input parameter signals that the system is 
always critical and one speaks of the system as being characterized by a "line 
of critical points." 

We already saw that U(l) gauge models with a lattice regulator confine in 

three Euclidean dimensions. All lattice gauge theories confine in two Euclidean 
dimensions. The demonstration is straightforward. Write down the partition 
function and gauge fix so that all the links in one direction (direction 2, in the 
following) are set to the identity. (It may be necessary to take open boundary 
conditions to do this.) The partition function becomes 



df7i(a;) exp(/3Tr[/i(x)C/|(a; + 2)). (71) 



Now perform a second set of gauge transformations Ui{x) = Ui{x)Ul{x + 2), 
etc. The partition function factorizes into a produce of one-plaquette actions 

Z=[j dU exp{(3TrU)]^'^ . (72) 
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The partition function has no singularities unless the single link integral does 
(this can happen at large-N, for example). It is easy to repeat the calculation 
for the Wilson loop and find an area law. 

Monte Carlo simulations show that four-dimensional U(l) gauge theory 
has a confining phase at strong coupling, apparently induced by monopole 
condensation as in three dimensions, and a phase transition at ~ 1 into a 
weak-coupling Coulomb phase. In that phase the relic monopole loops slightly 
renormalize the electric charge. The order of the transition was controversial 
for many years, but recent studies favor second order E3. If it were first order, 
the theory would have a continuum limit only in the Coulomb phase. If the 
transition is second order, one could define a continuum confined U(l) theory 
by approaching the phase transition from within the confined phase. 

Z(N) gauge theories in four dimensions have a phase diagram similar to 
spin models. There is a strong-coupling confinement phase, a magnetically- 
ordered Higgs phase at large (3, and at larger N, an intermediate coupling 
Coulomb phase. When there are only two phases, the deconfinement transition 
is first order. 

There are many cases where the lattice can be used to provide insight into 
a physical system, without the need for simulation. An example of such a 
situation is the gauge-Higgs system, with a (schematic) Lagrange density 

L^Lgauge + {D^,^f -V{^). (73) 

It is often useful to convert to radial and angular variables: $ = p0 with 
= 1. Then if we take V{^) = A($2 _ f^f = A(p2 _ pf^ we can freeze 
out the radial mode by taking the A — s- oo limit, so (p) = /. Replacing the 
gradient by a lattice difference, one can then write the action in terms of two 
coupling constants as 

L = l3gTr{l - UUUU) - MlU^{x)d^{x + fi). (74) 

The two coupling constants describe a plane of possible theories. One studies 
the continuum behavior of these theories on the lattice by finding all the critical 
points/lines in the phase diagram. One recovers continuum physics by tuning 
couplings to approach those critical points. 

One can usually figure out the behavior of the theory on the boundaries of 
the coupling constant plane without using Monte Carlo, and use Monte Carlo 
only to explore the interior of the phase diagram. This was a very common 
game in the days before the lattice community became "professional." Let's 
see how it works for a gauge-Higgs model. 

The first limiting case is (3h — 0. This is a pure gauge theory, which one 
has to analyze on its own. 
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Now Pg = 0. This is a trivial limit, the strong-coupling confined phase 
limit. One can see this in two ways. The first way is to realize that there is no 
gauge interaction, so all the link variables decouple and C = (p'^Ucj) describes 
a set of random interactions. All {4'i<tij) correlators vanish away from zero 
separation. Another way to see this limit is to gauge fix to "unitary gauge" by 
rotating all the 0's to some constant value 0o. Now 

Z = jY{dU^{x)e,qy{l3hY,4>lU^,{x)M (75) 

and again the partition function factorizes in to a product of one- link integrals. 

At /3g = oo we gauge-fix to an axial gauge by rotating all the [/'s to the 
identity. (Only U = 1 configurations are important in the functional integral 
in this limit.) Now — £ = Pi^cj)^ {x)(j){x + fi) which is a pure spin model. Again, 
one must analyze it on its own. 

Finally, we have /3h ^ oo. Use unitary gauge again. In this limit (pQUf^{x)(j>() 
is forced to its maximum value. If it happens the (j) transforms under a faithful 
representation of the gauge group G, only t7 = 1 is allowed. If the represen- 
tation of (f) is not faithful, then there is a subgroup H of G in which U ~ 1. 
Along this line, we then have a gauge theory with gauge group H (z G. pi the 
first case the Higgs and confinement phases are likely to be connected 113 (as 
shown in Fig. ^). In the second case, the gauge symmetry on the top line is 
typically a discrete symmetry, there is typically a first order phase transition 
(which extends into the phase diagram) and one must resort to Monte Carlo 
to see if the phases are connected or separated. Schematic phase diagrams for 
Z2 gauge- Z2 spin and gauge- Z3 spin models Ej are shown in Figs. ^ and 
b. 

Finally, a few words about fermions (see the next section for details). It is 
easy to study vectorlike theories on the lattice, and easiest to study multiples 
of 2 or 4 degenerate fiavors of fermions. It is reasonably easy to study the chiral 
properties of model theories, although taking the chiral limit in a computer 
simulation is hard. Common games to play are to measure (ipip) in a simulation 
at finite fermion mass and to extrapolate to zero mass. With Kogut-Susskind 
fermions, this operator is only multiplicatively renormalized, and so one can see 
if it vanishes at zero quark mass or extrapolates to some finite limit. One can 
also compute the pion mass and ask whether to J ~ iriq. Only fermion actions 
which are bilinear, S = ipMip for any matrix M, can be simulated on the 
computer, but if one has a four-fermion interaction, the Hubbard-Stratonovich 
transformation 

exp-(-0V')^= d(t)exp-{'ip->p(j) + 4>'^) (76) 
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Figure 4: Phase diagrams for (a) Z2 gauge-Z2 spin and and (b) Ze gauge-Zs spin models. 

allows four-fermion interaction models to be converted into Yukawa models 
(said differently, four-fermion interaction models are equivalent to Yukawa 
modelsllS), which the computer (and maybe the physicist) knows how to solve. 



5 Standard Methods for Spectroscopy 

5.1 Relativistic Fermions on the Lattice 

Defining fermions on the lattice presents a new problem: doubling. The naive 
procedure of discretizing the continuum fermion action results in a lattice 
model with many more low energy modes than one originally anticipated. Let's 
illustrate this with free field theory. 

The free Euclidean fermion action in the continuum (in four dimensions) 

is 

S = j (fx['^j{x)-i^d^,i}{x) + mi,{x)i,{x)]. (77) 

One obtains the so-called naive lattice formulation by replacing the derivatives 
by symmetric differences: we explicitly introduce the lattice spacing a in the 
denominator and write 

Sr''^ = ^«^(V'«+M - ^n-^) + mY,Mn. (78) 
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The propagator is: 



G{p) = {i-fu, sinp^a + ma) = — '— — (79) 

We identify the physical spectrum through the poles in the propagator, at 
Po = iE: 

sinh^ Ea = ^ sin^ pjU + vn?a?' (80) 



The lowest energy solutions are the expected ones at p = (0,0,0), E ~ ±m, 
but there are other degenerate ones, a.t p = (7r,0,0), (0,7r,0,), . . . (tt,??, tt). 
This is a model for sixteen light fermions, not one. 
(a) Wilson Fermions 

There axe two ways to deal with the doublers. The first way is to alter 
the dispersion relation so that it has only one low energy solution. The other 
solutions are forced to ~ 1/a and become very heavy as a is taken to zero. 
The simplest version of this solution (and almost the only one seen in the 
literature until recently) is due to Wilson: add a second-derivative-like term 

= ~^Y1 ^"(V'n+^ - 2V'n + V-n-^) (81) 

to 5"-"*"'=. The parameter r must lie between and 1; r = 1 is almost always 
used and "r = 1" is implied when one speaks of using "Wilson fermions." The 
propagator is 

^, , -n^sinp^a + ma-r^ (cosp^a-1) 
sm pi^a + {ma - r 2_^^{cosp^a - l)y 

It has one pair of poles at ~ (ztim, 0, 0, 0), plus other poles at p ~ r/a. In 
the continuum these states become infinitely massive and decouple (although 
decoupling is not trivial to prove). 

With Wilson fermions it is conventional not to use not the mass but the 
"hopping parameter" k = ^{ma + 4r)^^, and to rescale the fields tp ^/2Kip. 
The action for an interacting theory is conventionally written 

"5 = ^ i^nfpn - K^itpnir - 7^)[/^(n)V'n+^ + V'n(r + f^jUl^pn- fi) ■ (83) 
n nft 

Wilson fermions are closest to the continuum formulation- there is a four 
component spinor on every lattice site for every color and/or flavor of quark. 
Constructing currents and states is just like in the continuum. 



19 



However, the Wilson term explicitly breaks chiral symmetry. This has the 
consequence that the zero bare quark mass limit is not respected by interac- 
tions; the quark mass is additively renormalized. The value of Kc, the value of 
the (bare) hopping parameter at which the pion mass vanishes, is not known 
a priori before beginning a simulation; it must be computed. This is done in 
a simulation involving Wilson fermions by varying n and watching the pion 
mass extrapolate quadratically to zero as c^^ Kc — k {kc — k is proportional 
to the quark mass for small niq). 
(b) Staggered or Kogut-Susskind Fermions 

In this formulation one reduces the number of fermion flavors by using 
one component "staggered" fermion fields rather than four component Dirac 
spinors. The Dirac spinors are constructed by combining staggered fields on 
different lattice sites. 

Let us call the staggered field on site n Xn and define the matrix field 

^afe(n) = ^ir+%kXn+b (84) 

6 

where 

and b runs over the sites of a hypercube b = (0, 0, 0, 0), (1,0,0,0), . . . (1,1,1,1). 
•ipak is defined analogously, 

i'ak{n) = Y,ir+'tkXn+b- (86) 

b 

One interprets the a indices as spin (= 1 to 4) and k as flavor (= 1 to 4). 
The free action (with niq = 0) is 

^= ^nlMn+n - i'n-^i) (87) 

= ^^XnVn,n{Xn-n - Xn+n) (88) 

where ?7„^^ = (^^i^ni+n^+ns+rn _ Gauge fields can be added in the standard 
way, and explicit mass terms give an extra niq ^ XX term. 

Staggered fermions preserve an explicit chiral symmetry as niq — > even 
for finite lattice spacing, as long as all four flavors are degenerate. They are 
preferred over Wilson fermions in situations in which the chiral properties 
of the fermions dominate the dynamics-for example, in studying the chiral 
restoration/deconfinement transition at high temperature. They also present 
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a computationally less intense situation from the point of view of numerics than 
Wilson fermions, for the trivial reason that there are less variables. However, 
flavor symmetry and translational symmetry are all mixed together. Construc- 
tion of meson and haryon states (especially the A) is more complicated than 
for Wilson fermionsE3. 



5.2 Enter the Computer 

A "generic" Monte Carlo simulation in QCD breaks up naturally into two 
parts. In the "configuration generation" phase one constructs an ensemble 
of states with the appropriate Boltzmann weighting: we compute observables 
simply by averaging N measurements using the field variables 0^'^ appropriate 
to the sample 

1 ^ 

(r)^f ^-^r[</)«]. (89) 

i=l 

As the number of measurements N becomes large the quantity T will become 
|a^Gaussian distribution about a mean value. Its standard deviation is roughly 

N 

2 ' ' 



^F = ^(^Eir['/'^*^]l'-r2)- (90) 



=1 



The idea of essentially all simulation algorithms is that one constructs a new 
configuration of field variables from an old one. One begins with some simple 
field configuration and monitor observables while the algorithm steps along. 
After some number of steps, the value of observables will appear to become 
independent of the starting configuration. At that point the system is said to 
be "in equilibrium" and Eq. ^ can be used to make measurements. 
i— I The simplest method for generating configurations is called the Metropolis 
113 algorithm. It works as follows: From the old configuration {(/)} with action 
(3S, transform the variables (in some reversible way) to a new trial config- 
uration {(/)}' and compute the new action /3S'. Then, if S" < make the 
change and update all the variables; if not, make the change with probability 
eM-PiS'-S)). 

Why does it work? In equilibrium, the rate at which configurations i turn 
into configurations j is the same as the rate for the back reaction j ^ i. The 
rate of change is (number of configurations) x (probability of change). Assume 
for the sake of the argument that Si < Sj. Then the rate i ^ j is NiP{i j) 
with P{i — > j) = exp(— /3(S'j — Si) and the rate j — > i is NjP{j — > i) with 
P{j ^i) = l. Thus N,/Nj = exp(-/3(S', - Sj)). 
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If you have any interest at all in the techniques I am describing, you 
should write a little Monte Carlo program simulating the two-dimensional Ising 
model. Incidentally, the favorite modern method for pure gauge models is 
overrelaxationL3. 

One complication for QCD which spin models don't have is fermions. The 
fermion path integral is not a number and a computer can't simulate fermions 
directly. However, one can formally integrate out the fermion fields. For n/ 
degenerate flavors of staggered fermions 

r _ _ 

Z= [dU] [d^P] [d^:] eM-l3S{U) - ^ ^pMifj) (91) 

i=i 

[dU](detM)"'f^^exp{-pS{U)). (92) 

(One can make the determinant positive-definite by writing it as det(M'''M)"^ /^.) 
The determinant introduces a nonlocal interaction among the t/'s: 

Z = J [dU] exp{-f3S{U) - ^Trln(Af^Af)). (93) 

All large scale dynamical fermion simulations today generate configura- 
tions using some variation of the microcanonical ensemble. That is, they in- 
troduce momentum variables P conjugate to the J7's and integrate Hamilton's 
equations through a simulation time t 

if = iPU (94) 

The integration is done numerically by introducing a timestep At. The mo- 
menta are repeatedly refreshed by bringing them in contact with a heatj_bath 
and the method is thus called Refreshed or Hybrid Molecular DynamicsEj. 

For special values of rt/ (multiples of 2 for Wilson fermions or of 4 for 
staggered fermions) the equations of motion can be derived from a local Hamil- 
tonian and in that case At systematics in the integration can be removed by 
an extra Metropolis accept/reject step. This method is called Hybrid Monte 
Carlo0. 

The reason for the use of these small timestep algorithms is that for any 
change in any of the ?7's, (M'l'M)^^ must be recomputed. When Eqn. ^ is 
integrated all of the ?7's in the lattice are updated simultaneously, and only 
one matrix inversion is needed per change of all the bosonic variables. 
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The major computational problem dynamical fermion simulations face is 
inverting the fermion matrix M. It has eigenvalues with a very large range- 
from 2ti down to niqa- and in the physically interesting limit of small ruq 
the matrix becomes ill-conditioned. At present it is necessary to compute at 
unphysically heavy values of the quark mass and to extrapolate to rUq = 0. 
The standard inveraipn technique today is one of the variants of the conjugate 
gradient algorithm l3. 

5.3 Spectroscopy Calculations 

"In a valley something like a race took place. A little crowd watched bunches of 
cars, each consisting of two 'ups' and a 'down' one, starting in regular intervals 
and disappearing in about the same direction. 'It is the measurement of the 
proton mass,' commented Mr. Strange, 'thej2.have done it for ages. A very 
dull job, I am glad I am not in the game.' " 

Masses are computed in lattice simulations from the asymptotic behavior 
of Euclidean-time correlation functions. A typical (diagonal) correlator can be 
written as 

C{t) = (0|OWO(0)|0). (96) 

Making the replacement 

0{t) = e^*Oe-^* (97) 
and inserting a complete set of energy eigenstates, Eq. (3.1) becomes 

C(0=5]|(0|O|n)pe-^"*. (98) 

n 

At large separation the correlation function is approximately 

C(i) ^ |(0|O|l)pe-^i* (99) 

where Ei is the energy of the lightest state which the operator O can create 
from the vacuum. If the operator does not couple to the vacuum, then in the 
limit of large t one hopes to to find the mass Ei by measuring the leading 
exponential falloff of the correlation function, and most lattice simulations 
begin with that measurement. If the operator O has poor overlap with the 
lightest state, a reliable value for the mass can be extracted only at a large 
time t. In some cases that state is the vacuum itself, in which Ei = 0. Then 
one looks for the next higher state-a signal which disappears into the constant 
background. This makes the actual calculation of the energy more difhcult. 

This is the basic way hadronic masses are found in lattice gauge theory. 
The many calculations differ in important specific details of choosing the op- 
erators 0{t). 
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5-4 Recent Results 



Today's supercomputer QCD simulations range from 16"^ x 32 to 32"^ x 100 
points and run from hundreds (quenched) to thousands (fuU QCD) of hours 
on the fastest supercomputers in the world. 

Results are presented in four common ways. Often one sees a plot of some 
bare parameter vs. another bare parameter. This is not very useful if one wants 
to see continuum physics, but it is how we always begin. Next, one can plot a 
dimensionless ratio as a function of the lattice spacing. These plots represent 
quantities like Eqn. Both axes can show mass ratios. Examples of such 
plots are the so-called Edinburgh plot {mN/rrip vs. m-n/mp) and the Rome 
plot {misi/mp vs. {mT^/mp)'^). These plots can answer continuum questions 
(how does the nucleon mass change if the quark mass is changed?) or can be 
used to show (or hide) scaling violations. Plots of one quantity in MeV vs. 
another quantity in MeV are typically rather heavily processed after the data 
comes off the computer. 

Let's look at some examples of spectroscopy, done in the "standard way," 
with the plaquette gauge action and Wilson or staggered quarks. I will restrict 
the discussion to quenched simulations because only there are the statistical 
errors small enough to be interesting to a non-lattice audience. Most dynamical 
fermion simulations are unfortunately like Dr. Johnson's dog. 

Fig. H shows a plot of the rho mass as a function of the size of the simula- 
tion, for several values of the quark mass (or m^/mp ratio in the simulatioiL) 
and lattice spacing (/3 = 6.0 is a ~ 0.1 fm and (3 = 5.7 is about twice that)Ej. 
This picture shows that if the box has a diameter bigger than about 2 fm, the 
rho mass is little affected, but if the box is made smaller, the rho is "squeezed" 
and its mass rises. 

Next we look at an Edinburgh plot. Fig. The different plotting 

symbols correspond to different bare couplings or (equivalently) different lattice 
spacings. This plot shows large scaling violations: mass ratios from different 
lattice spacings do not lie on top of each other. We can expose the level of 
scaling violations by taking "sections" through the plot and plot tojv/^p at 
fixed values of the quark mass (fixed m-^jmp)^ vs. lattice spacing, in Fig. ^ 

Now for some examples of scaling tests in the chiral limit. (Extrapolating 
to the chiral limit is a whole can of worms on its own, but for now let's assume 
we can do it.) Fig. || shows the nucleon to rho niass ratio (at diiral limit) 
vs. lattice spacing (in units of l/wp) for staggeredEj and WilsonES fermions. 
The "analytic" result is from strong coupling. The two curves are quadratic 
extrapolations to zero lattice spacing using different sets of points from the 
staggered data set. The burst is from a linear extrapolation to the Wilson 
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Figure 6: An Edinburgh plot for staggered fermions, from tiie MILC collaboration. 
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data. The reason I show this figure is that one would like to know if the 
continuum limit of quenched spectroscopy "predicts" the real- world N/p mass 
ratio of 1.22 or not. The answer (unfortunately) depends on how the reader 
chooses to extrapolate. 

Another tested is the ratio of the rho mass to the square root of the string 
tension, Fig. ^. Here the diamonds are staggered data and the crosses from 
the Wilson action. Scaling violations are large but the eye extrapolates to 
something close to data (the burst). 

Finally, despite Mr. Strange, very few authors have attempted to extrap- 
olate to infinite volume, zero lattice spacing, and to physical quark masses. 



including the strange quark. One group which did, Butler et al. c3, produced 



Fig. 10, The squares are lattice data, the octagons are the real world. They 
look quite similar within errors. Unfortunately, to produce this picture, they 
had to build their own computer. 



6 Bare Actions, Scaling, and Improved Actions 

(Disclaimer: It happens that at present there are a number of approaches 
to "improvement." At present, the words which are used to describe them 
differ greatly among the different groups, and the physics motivations are quite 
different. For example, in one approach, ideas based on perturbation theory are 
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Figure 8; Nucleon to rho mass ratio (at chiral limit) vs. lattice spacing (in units of l/mp). 
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Figure 9: Scaling test for the rho mass in terms of the string tension, with data points labeled 

as in Fig. 8. 
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Figure 10: Quenched approximation mass ratios from Ref. 25. 



regarded as absolutely crucial while in another approach perturbation theory 
is said to be completely irrelevant. The particular problems whose solution 
would be regarded as a success are different for different groups. This makes 
life exciting for the researcher but dangerous for the reviewer. You should 
know at the start that I work in this area, and that I am as prejudiced as 
anyone else in the field.) 

The slow approach to scaling presents a practical problem for QCD sim- 
ulations, since it means that one needs to work at small lattice spacing. This 
is expensive. The cost of a Monte Carlo simulation in a box of physical size L 
with lattice spacing a and quark mass ruq scales roughly as 

i-n-y-'i—r-' im 



a a m, 



q 



where the 4 is just the number of sites, the 1-2 is the cost of "critical slowing 
down" the extent to which successive configurations are correlated, and the 
2-3 is the cost of inverting the fermion propagator, plus critical slowing down 
from the nearly massless pions. The problem is that one needs a big computer 
to do anything. 

However, all the simulations I described in the last lecture were done with 
a particular choice of lattice action: the plaquette gauge action, and either 
Wilson or staggered quarks. While those actions are the simplest ones to 
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program, they are just particular arbitrary choices of bare actions. Can one 
invent a better lattice discretization, which has smaller scaling violations? 

To approach this problem, let's think about the connection between scaling 
and the properties of some arbitrary bare action, which we assume is defined 
with some UV cutoff a (which does not have to be a lattice cutoff, in principle). 
The action is characterized by an infinite number of coupling constants, {g}. 
When the g^s have any arbitrary value, the typical scale for all physics will 
be the order of the cutoff: m ~ f/a, correlation length ^ ~ a. There will be 
strong cutoff effects. 

_ The best way to think about scaling is through the renormalization group 
EZl. Take the action with cutoff a and integrate out degrees of freedom to 
construct a new effective action with a new cutoff a' > a and fewer degrees 
of freedom. The physics at distance scales r > a is unaffected by the coarse- 
graining (assuming it is done exactly.) We can think of the effective actions as 
being similar to the original action, but with altered couplings. We can repeat 
this coarse-graining and generate new actions with new cutoffs. As we do, the 
coupling constants "flow:" 

S(a, cj) S{2a, c^-) ^ S{ia, c^') ^ ... (101) 

If under repeated blockings the system flows to a fixed point 

5(a„, c]) ^ S{a^+i,c]+^ = c]) (102) 

then observables are independent of the cutoff a and in particular the correla- 
tion length ^ must either be zero or infinite. 

This can only happen if the original c's belong to a particular restricted 
set, called the "critical surface." It is easy to see that physics on the critical 
surface is universal: at long distances the effective theory is the action at the 
fixed point, to which all the couplings have flowed, regardless of their original 
bare values. 

But ^ = cxD is not f large. Imagine tuning bare parameters close to the 
critical surface, but not on it. The system will flow towards the fixed point, 
then away from it. The flow lines in coupling constant space will asymp- 
totically approach a particular trajectory, called the renormalized trajectory, 
which connects (at ^ = 00) with the fixed point. Along the renormalized tra- 
jectory, ^ is finite. However, because it is connected to the fixed point, it shares 
the scaling properties of the fixed point-in particular, the complete absence of 
cutoff effects in the spectrum and in Green's functions. (To see this remarkable 
result, imagine doing QCD spectrum calculations with the original bare action 
with a cutoff equal to the Planck mass-or beyond (!) and then coarse graining. 
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Figure 11: A schematic picture of renormalization group flows along a one-dimensional 
critical surface, with the associated renormalized trajectory, and superimposed contours of 

constant correlation length. 



Now exchange the order of the two procedures. If this can be done without 
making any approximations the answer should be the same.) 

A Colorado analogy is useful: think of the critical surface as the top of 
a high mountain ridge. The fixed point is a saddle point on the ridge. A 
stone released on the ridge will roll to the saddle and come to rest. If it is not 
released exactly on the ridge, it will roll near to the saddle, then go down the 
gully leading away from it. For a cartoon, see Fig. |ll|. 

So the ultimate goal of "improvement programs" is to find a true per- 
fect action, without cutoff effects, along the renormalized trajectory of some 
renormalization group transformation. 

The bare action for an asymptotically free theory like QCD can be char- 
acterized by one weakly relevant operator, F^^,, and its associated coupling 
5^, plus many irrelevant operators. In lattice language, our bare actions are 
described by one overall factor of /3 = "^N/ and arbitrary weights of various 
closed loops, 

PS -^11^=0,. (103) 

Asymptotic freedom is equivalent to the statement that the critical surface of 
any renormalization group transformation is at = 0. The location of a fixed 
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point involves some relation among the Cj's. 

As you might imagine, finding a renormalized trajectory is hard, and no- 
body has done it in a convincing way yet for QCD. What are people doing in 
the meantime? 

6.1 Improvement based on naive dimensional analysis 

The idea here is that since the critical surface is at ^ we can use the naive 
canonical dimensionality of operators to guide us in our choice of improvement. 
If we perform a naive Taylor expansion of a lattice operator like the plaquette, 
we find that it can be written as 

1 - ^Re TiUpiaq = roTiF^, + [n "^^D^F^^.D^F^, + 

^2E,,..Tri?^^.-^M^- + 
^3E^..Tri^M^^,i?,F,,] + 

-fO(a'*) (104) 

The expansion coefficients have a power series expansion in the coupling, rj = 
Aj + g^Bj + . . . and the expectation value of any operator T will have an 
expansion 

(r(a)) = (r(0)) + 0(a) + 0{g\) + ... (105) 

Other loops have a similar expansion, with different coefficients. Now 
the idea is to take a minimal subset of loops and systematically remove the 
a" terms order by order for physical observables by taking the right linear 
combination of loops. 

5 = 5]c,0, (106) 

3 

with 

c,=c°+g^ + ... (107) 

"Tree level" improvement removes pure power-law corrections through some 
order in a". One can also consider quantum corrections, and at each order 
in a" remoizK-oIIL terms, too. This method was developed by Symanzik and 
co-workersc3E3 l2I ten years ago. 

The ipDst commonly used "improved" fermion action is the "Sheikholeslami- 
Wohlert"E2l or "clover" action, an order improved Wilson action. The orig- 
inal Wilson action has 0{a) errors in its vertices, Sw = Sc + 0{a). This is 
corrected by making a field redefinition 

i;{x)^ib'{x)^ i;{x) + fj^D^i; (108) 
^Pix)^^'{x)= ^Pix) + f-f^^D^ (109) 
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Figure 12: The "clover term" 



and the net result is an action with an extra lattice anomalous magnetic mo- 
ment term, 

Ssw - '^-^^ix)(Tf,^Ff,^^{x) (110) 
It is called the "clover" action because the lattice version of F^^i, is the sum of 



paths shown in Fig. 12 



Studies performed at the time showed that this program did not improve 
scaling for the pure gauge theory (in the sense that the cost of simulating 
the more complicated action was greater than the savings from using a larger 
lattice spacing.) The whole ppagram was re-awakened in the last few last years 
by Lepage and collaborators 123 and variations of this program give the most 
widely used "improved" lattice actions. 



6.2 Nonperturbative determination of coefficients 

Although I am breaking chronological acjder, the simplest approach to Symanzik 
improvement is the newest. The ideaH is to force the lattice to obey various 
desirable identities to some order in a, by tuning parameters until the identities 
are satisfied by the simulations. Then calculate other things. One example is 
the PCAC relation 

9^^;^ = 2m?, + 0(a), (111) 
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where the axial and pseudoscalar currents are just 



Al{x)=ij{x)^^^^-T-i>{x) (112) 

and 

P-{x) = V;(x)75iT'^7/;(x) (113) 
(r" is an isospin index.) The PCAC relation for the quark mass is 

0{a). (114) 



2 (P'-C) 



Now the idea is to take some Symanzik-improved action, with the expansion 
coefficients allowed to vary, and perform simulations in a little box with some 
particular choice of boundary conditions for the fields. Tune the parameters, 
which include the «/4c5i4'Ocr^i/P^t^ and redefined currents 

A^ = ZA[{l + bAam^)A''^ + CAad^P'' (115) 
P" = Zp{l + bpamg)P'', (116) 



until the quark mass, defined in Eq. 114, is independent of location in the box, 
or of the boundary conditions. Figs. |13| and ^ illustrate what can be done 
with this tuning procedure. It is still too soon for definitive tests of scaling 
with this procedure. 



6.3 Improving perturbation theory 

The older version of Symanzik improvement uses lattice perturbation theory 
to compute the coefficients of the operators in the action. Let's make a di- 
gression into lattice perturbation theory lJ. It has three major uses. First, 
we need to relate lattice quantities (like matrix elements) to continuum ones: 
0™"*(/x) = Z(/xa, (7(a))0'°**(a). This happens because the renormalization of 
an operator is slightly different in the two schemes. In perturbation theory 
Z has an expansion in powers of g^. Second, we can use perturbation theory 
to understand and check numerical calculations when the lattice couplings are 
very small. Finally, on£_can use perturbative ideas to motivate nonperturbative 
improvement schemes l3. 

Perturbation theory for lattice actions is just like any other kind of per- 
turbation theory (only much messier). One expands the Lagrangian into a 
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Figure 13: Values of the quark mass as computed from the axial and pseudoscalar cur- 
rents, using the Wilson action. The open and full symbols correspond to different boundary 

conditions on the gauge fields. 




Figure 14: Same as previous figure, but now with improved action and operators. 
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quadratic term and interaction terms and constructs the propagator from the 
quadratic terms: 

C= A^{x)pf,^{x-y)A^{y) + gA^ + ... (117) 
Lo+Li. (118) 

For example, the gluon propagator in Feynman gauge for the Wilson action is 

Ep(l -cos(g^a)) 

To do perturbation theory for any system (not just the lattice) one has to 
do three things: one has to fix the renormalization scheme (RS) (define a 
coupling), specify the scale at which the coupling is defined, and determine a 
numerical value for the coupling at that scale. All of these choices are arbitrary, 
and any perturbative calculation is intrinsically ambiguous. 

Any object which has a perturbative expansion can be written 

0{Q) = CO + ci(Q//i, RS)as{ti, RS) + c2(Q/m, RS)as{ti, RSf + ... (120) 

In perturbative calculations we truncate the series after a fixed number of 
terms and implicitly assume that's all there is. The coefficients Ci{Q/fj,,RS) 
and the coupling as{p, RS) depend on the renormajliaation scheme and choice 
of scale fi. The guiding rule of perturbation theory Ell is "For a good choice of 
expansion the uncalculated higher order terms should be small." A bad choice 
has big coefficients. 

There are many ways to define a coupling: The most obvious is the bare 
coupling; as we will see shortly, it is a poor expansion parameter. Another 
possibility is to define the coupling from some physical observable. One popular 
choice is to use the heavy quark potential at very high momentum transfer to 
define 

Viq)^4nC,^. (121) 

There are also several possibilities for picking a scale: One can use the bare 
coupling, then p — 1/a the lattice spacing. One can guess the scale or or 
play games just like in the continuum. One game is the Lepage-Mackenzie q* 
prescription: take 

".(-Z*) / Ae(<z) - / d^qmo^siq). (122) 



To find q* , write Ois{q) = ots{p) + &ln(g //i )as(/i) + . . ., and similarly for 



as{q*), insert these expressions into Eq. 122 and compare the as(/^)^ terms 
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to get 

Hq*) = I d\Hq)(/ J d\^. (123) 

This is the lattice analog of the Brodsky-Lepage-Mackenzie0 prescription in 
continuum PT. 

Finally one must determine the coupling: If one uses the bare lattice cou- 
pling it is already known. Otherwise, one can compute it in terms of the bare 
coupling: 

ajjg{s/a) = ao + (5.88-1.751ns)ag+(43.41-21.891ns+3.061n2s)ai] + . . . (124) 

Or one can determine it from something one measures on the lattice, which 
has a perturbative expansion. For example 

1 47r 
- \n{-TiUpiaq) = —ap{iAl/a){l - l.lSSap) (125) 

(to this order, ap — ay). Does "improved perturbation theory" actually 
improve perturbative jCalculations? In many cases, yes: some examples are 
shown in Fig. [l^ fromtJ: On the upper left we see a calculation of the average 
link in Landau gauge, from simulations (octagons) and then from lowest-order 
perturbative calculations using the bare coupling (crosses) and ay and a^/g 
(diamonds and squares). In the upper right panel we see how the lattice 
prediction of an observable involving the 2 by 2 Wilson loop depends on the 
choice of momentum q* / a (at (3 = 6.2, a rather weak value of the bare coupling) 
in the running coupling constant. The burst is the value of the prescription of 



Eq. 123 . In the lower panel are perturbative predictions the same observables 
as a function of lattice coupling. These pictures illustrate that perturbation 
theory in terms of the bare coupling does not work well, but that using other 
definitions for couplings, one can get much better agreement with the lattice 
"data" . 

Straight perturbative expansions by themselves for the commonly-used 
lattice actions are typically not very convergent. The culprit is the presence of 
?7^'s in the action. One might think that for weak coupling, one could expand 

{l)U^ = 4j[l + igaA + (126) 

and ignore the . . ., but the higher order term ■4)^g'^a? A^ij] generates the "tad- 
pole graph" of Fig. |l^. The UV divergence in the gluon loop ~ l/o? cancels 
the a? in the vertex. The same thing happens at higher acder, and the tad- 
poles add up tO|-ap effective aP ^ Cng^^ contribution. ParisicS and later Lepage 
and Mackenzie ca suggested a heuristic way to deal with this problem: replace 



36 



0.3 



0.2 



^ 0.1 



0.0 



5.5 



ite Carlo o 

av(q') * 
gMs(q*) ° 

a o 

D 



6.0 



6.5 



0.3 



0.2 — 



0.1 — 



0.0 



(b) 



2 
q/a 



0.4 



0.2 



0.0 









□ 






: X X 


X 




- «lat X 




- av(q*) 






; «Ms(q') n 






Monte Carlo o 







5.6 



6.2 6.4 



5.8 6.0 

Figure 15: Examples of "improved perturbation theory" 
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Figure 16: The "tadpole diagram" 



uq{1 + igaA) where uq, the "mean field term" or "tadpole improvement 
term" is introduced phenomenologically to sum the loops. Then one rewrites 
the action as 

^=iETrC/--J-ETrU (127) 
9 y/at 

where = gfat/'^o ^'^^ expansion parameter. Is Uq = (TrC/p/ag/S)? This 

choice is often used; it is by no means unique. 

A "standard action" (for this year, anyway) is the "tadpole-improved 
Liischer- Weisz action," composed of a 1 by 1, 1 by 2, and "twisted" loop 

(+x,+y,+z,-x,-y,-z), 

pS = -/?[Tr(l X 1) - —^(1 + 0.48a,)Tr(l x 2) - ^0.33a,Trf/t.„] (128) 

with Uq = (Trl/p;a<j/3)^/'* and 3.068as = —ln{TTUpiaq/3) determined self- 
consistently in the simulation. 

None of this section seems to have anything to do with critical surfaces or 
renormalized trajectories. What is going on? I believe that the connection is 
that people following this path to improvement are (implicitly) trying to find 
a trajectory in coupling constant space for which physical observables have 
no cutoff effects through some order in a" and/ or (g ^)'"Q". The trajectory 
is described by the variation of the Cj's in Eqn. |106| with respect to a or to 
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bare g^. If perturbation theory is reliable, one could compute them using the 
expansion of Eqn. 107. This is not a renormalized trajectory, but it might be 



"good enough" in the engineering sense. 



6.4 Fixed Point Actions 

A more direct attack on the renormalized trajectory begins by finding a fixed 
point action. Imagine having a set of field variables {0} defined with a cutoff 
a. Introduce some coarse-grained variables {$} defined with respect to a new 
cutoff a', and integrate out the fine-grained variables to produce a new action 

e-/3S'(*) = j dcl,e-^(n^'^)+SM) (129) 
where /3(T(<I>, 0) is the blocking kernel which functionally relates the coarse and 



fine variables. Integrating Eq. 12£ is-.usually horribly complicated. However, 
P. Hasenfratz and F. Niederniayer l3 noticed an amazing simplification for 
asymptotically free theories: Their critical surface is at /? = 00 and in that 



limit Eq. 129 becomes a steepest-descent relation 

^'($) = min((r($, (j)) + 5(</))) (130) 



which can be used to find the fixed point action 



Sfp{^) = min((r($, 0) + Sppm. (131) 

The program has been successfully carri£jd out for d ~ 2 sigma models El and 
for four-dimensional pure gauge theoriesEj. These actions have two noteworthy 
properties: First, not only are they classically perfect actions (they have no 
a" scaling violations for any n), but they are also one-loop quantum perfect: 
that is, as one moves out the renormalized trajectory, 

\srt{9^) = \{Sfp + 0{g^)). (132) 

g g 

Physically this happens because the original action has no irrelevant operators, 
and they are only generated through loop graphs. Thus these actions are an 
extreme realization of the Symanzik program. Second, because these actions 
are at the fixed point, they have scale invariant classical solutions. This fact 
can be used to define a topological charge, operator on the lattice in a way 
which is consistent with the lattice actionEa. 
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These actions are "engineered" in the foUowing way: one picks a favorite 



blocking kernel, which has some free parameters, and solves Eq. 131, usually 
approximately at first. Then one tunes the parameters in the kernel to optimize 
the action for locality, and perhaps refines the solution. Now the action is used 
in simulations at finite correlation length (i.e. do simulations with a Boltzman 
factor exp{—pSFp)- Because of Eq. 132, one believes that the FP action will 
be a good approximation to the perfect action on the RT; of course, only a 
numerical test can tell. As we will see in the next section, these actions perform 
very well. At this point in time no nonperturbative FP actioiL which includes 
fermions has been tested, but most of the formalism is there EJ. 



6.5 Examples of "Improved" Spectroscopy 

I would like to show some examples of the various versions of "improvement" , 
and remind you of the pictures at the end of the last chapter to contrast results 
from standard actions. 

Fig. |l^ shows a plot of the string tension measured in systems of constant 
physical size (measured in units of 1/Tc, the critical temperature for dccon- 
finement), for SU(3) pure gauge theory. In the quenched approximation, with 

~ 440 MeV, Tc = 275 MeV and l/Tc = 0.7 fm. Simulations with the 
standard Wilson action are crosses, while the squares show FP action-cesults 
and the octagons from the tadpole- improved Liischer-Weisz action e3. The 
figure illustrates that it is hard to quantify improvement. There are no mea- 
surements with the Wilson action at small lattice spacing of of precisely the 
same observables that the "improvement people" measured. The best one can 
do is to take similar measurements (the diamonds) and attempt to compute 
the a = prediction for the observable we measured (the fancy cross at a = 0). 
This attempt lies on a straight line with the FP action data, hinting strongly 
that the FP action is indeed scaling. The FP action seems to have gained 
about a factor of three to four in lattice spacing , or a gain of (3 — 4)^ com- 
pared to the plaquette action, according to Eq. |lO0| , at a cost of a factor of 
7 per site because it is more complicated to program. The tadpole-improved 
Liischer-Weisz action data lie lower than the FP action data and do not scale 
as well. As a ^ the two actions should yield the same result; that is just 
universality at work. However, there is no guarantee that the approach to the 
continuum is monotonic. 

Fig. ^ shows the heavy quark- antiquark potential in SU(2) gauge theory, 
where V{r) and r are measured in the appropriate units of Tc, the critical 
temperature for deconfinement. The Wilson action is on the left and a FP 
action is on the right. The vertical displacements of the potentials are just 
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Figure 17: The (square root of) the string tension in lattices of constant physical size 
L = 2/Tc, but different lattice spacings (in units of l/Tc). 



there to separate them. Notice the large violations of rotational symmetry 
in the Wilson action data when the lattice spacing is a = l/2Tc which are 
considerably improved in the FP action results. 

Next we consider nonrelativi^tic QCD. A comparison of the quenched char- 
monium spectrum from Rcf. L; 



19. When the tadpole- 



is shown in Fig. 

improved L-W action is used to generate gauge configurations, the scaling 
window is pushed out to a ~ 0.4 fm for these observables. 

Now we turn to tests of quenched QCD for light quarks. The two actions 
which have been most extensively tested are the S-W action, with and without 
tadpole improvement, and an axition called the D234(2/3) action, a higher- 
order variant of the^-W actionC3. J'igs. ^ and ^ are the analogs of Figs. 
^ and H Diamonds c3 and plussesc3 are S-W actions, ordinary and tadpole- 
improved, squares are the D234(2/3) action. They appear to have about half 
the scaling violations as the standard actions but they don't remove all scaling 
violations. It's a bit hard to quantify the extent of improvement from these 
pictures because a chiral extrapolation is hidden in them. However, one can 
take one of the "sections" of Fig. |^ and overlay the new data on it. Fig. It 
looks like one can double the lattice spacing for an equivalent amount of scale 
violation. However, the extrapolation in a is not altogether clear. Fig. ^ is 
the same data as Fig. p3, only plotted vs. a^, not a. All of the actions shown 
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Figure 18: The Heavy quark potential in SU(2) pure gauge theory measured in units of T^. 

(a) Wilson action (b) an FP action. 
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Figure 19: S, P, and D states of charmonium computed on lattices with: a = 0.40 fm 
(improved action, Ppiaq = 6-8); a = 0.33 fm (improved action, Ppi^q = 7.1); a = 0.24 fm 
(improved action, (ipiaq = 7.4); and a = 0.17 fm (Wilson action, /3 = 5.7, from [43]), from 
Ref. 32. The dashed lines indicate the true masses. 
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Figure 20: Nuclcon to rho mass ratio (at chiral limit) vs. lattice spacing (in units of l/nip). 



in these figures are supposed to have 0{a?) (or better) scaling violations. Do 
the data look any straighter in Fig. psl than Fig. 22 ? 



6.6 The bottom line 

At the cost of enormous effort, one can do fairly high precision simulations 
of QCD in the quenched approximation with standard actions. The actions I 
have shown you appear to reduce the amount of computation required for pure 
gauge simulations from supercomputers to very large work stations, probably a 
gain of a few hundreds. All of the li ght quark data I showed actually came from 



supercomputers. According to Eq. IOC , a factor of 2 in the lattice spacing gains 



a factor of 64 in speed. The cost of either of the two improved actions I showed 
is about a factor of 8-10 times the fiducial staggered simulation. Improvement 
methods for fermions are a few years less mature than ones for pure gauge 
theory, and so the next time you hear a talk about the lattice, things will have 
changed for the better (maybe). 



7 Case Studies 
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Figure 21: Rho mass scaling test with respect to the string tension. 
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Figure 22: miv/mp vs amp at fixed quark mass (fixed TriTr/mp). Interpolations of the S-W 
and D234(2/3) data were done by me. 



44 



L8 



1-6 



1.4 



1-2 



1.0 

1 2 

(am^(a = 0))" 

Figure 23: nifq jnip vs (arrip)^ at fixed quark mass (fixed m-n/mp). 



7.1 Glueballs 

To this audience, glueballs are just the closed strings of QCD. They are in- 
teresting because they are part of the QCD spectrum but not part of a naive 
quark model. What do we expect for a spectrum? Models like the bag model, 
which have bound "constituent gluons," would predict that the lightest state is 
a scalar, followed by tensor and pseudoscalar states, but nobody really knows. 
Here is a place where the lattice is the only game in town for a first-principles 
calculation. 

Since this is a TASI for stringy graduate students, the experimental sit- 
uation can be summarized briefly: it is a mess c3. There are two leading 
experimental candidates, the /j(1700), seen in radiative ip decays {ip jKK 
and ip ^rjrj (mostly at SLAC)) and the /o(1500), seen in pp T^^rj-q, from 
the Crystal Barrel detector at CERN. 

People have been trying to measure the masses of the lightest glueballs (the 
scalar and the tensor) using lattice simulations for many years. The problem 
has proven to be very hard, for several reasons. 

Recall how we measure a mass from a correlation function (Eqn. |9^) . The 
problem with the scalar glueball is that the operator O has nonzero vacuum 
expectation value, and the correlation function approaches a constant at large 
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lim C(i)-> |(0|O|p = 0)pexp(-TOt) + |(0|O|0)p. (133) 

t — >oo 

The statistical fluctuations on C{t) are given by Eq. |9^ and we find after a 
short calculation that 

~ . (134) 



N 

Thus the signal to noise ratio collapses at large t like vW exp(— mt) due to the 
constant term. 

A partial cure for this problem is a good trial wave function O. While in 
principle the plaquette itself could be used, it is so dominated by ultraviolet 
fluctuations that it does not produce a good signal. Instead, people invent 
"fat links" which average the gauge field over several lattice spacings, and 
then make interpolating fields which are closed loops of these fat links. The 
lattice glueball is a smoke ring. 

The second problem is that lattice actions can have phase transitions at 
strong or intermediate coupling, which_have nothing to do with the continuum 
limit, but mask continuum behavior c3. As an example of this, consider the 
gauge group SU{2), where a link matrix can be parameterized as C/ = 1 cos 9 + 
ia-ns'm 6, so TiU — 2 cos 9. Now consider a generalization of the Wilson action 
—S = f3TrU + 7(TrJ7)^. (this is a mixed fundamental-adjoint representation 
action). At 7 — > cxd TtU — > ±1 and the gauge symmetry is broken down to 
Z{2). But Z{2) gauge theories have a first order phase transition. First order 
transitions are stable under perturbations, and so the phase diagram of this 



theory, shown in Fig. 24 , has a line of first order transitions which terminate in 
a second order point. At the second order point some state with scalar quantum 
numbers becomes massless. However, now imagine that you are doing Monte 
Carlo along the 7 = line, that is, with the Wilson action. When you come 
near the critical point, any operator which couples to a scalar particle (like the 
one you are using to see the scalar glueball) will see the nearby transition and 
the lightest mass in the scalar channel will shrink. Once you are past the point 
of closest approach, the mass will rise again. Any scaling test which ignores 
the nearby singularity will lie to you. 

This scenario has been mapped out for SU{3), and the place of closest 
approach is at a Wilson coupling corresponding to a lattice spacing of 0.2 
fm or so, meaning that very small lattice spacings are needed before one can 
extrapolate to zero lattice spacing. A summary of the sitiiation is shown in 
Fig. Here the quantity tq is the "Sommer radius" l3, defined through 

the force, by r§i^(ro) = —1.65. In the physical world of three colors and four 
flavors, ro = 0.5 fm. 
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Figure 24: Phase transitions in the fundamental-adjoint pli 
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Figure 25; Glueball mass vs ro with the Wilson action, from a summary picture in Ref. 49. 




47 



Finally, other arguments suggest that a small lattice spacing or a good 
approximation to an action on an RT are needed to for glueballs; the physical 
diameter of the glueball, as inferred from the size of the best interpolating 
field, is small, about 0.5 fm. Shafer and ShuryakEd have argued that the small 
size is due to instanton effiects. Most lattice actions at large lattice spacing do 
bad things to instantonsc3. 

Two big simulations have carried calculations oLthe glueball mass close 
to the continuum limit: the UKQGD collaborations and a collaboration at 
IBM which built its own computer E3. (The latter group is the one with the 
press release last December announcing the discovery of the glueball.) Their 
predictions in MeV are different and they each favor a different experimental 
candidate for the scalar glueball (the one which is closer to their prediction, 
of course). It is a useful object lesson because both groups agree that their 
lattice numbers agree before extrapolation, but they extrapolate differently to 
a = 0. 

The UKQCD group sees that the ratio m{Q'^^)/^a can be well fitted 
with a form h + ca?a [a is the string tension) and a fit of this form to the 
lattice data of both groups gives m(0+"'")/-\/cr = 3.64 ± 0.15. To turn this into 
MeV we need a in MeV units. One way is to take mpj^a and extrapolate 
that to a = using h + ca^a. Averaging and putting 770 MeV for mp one 
gets \Ju = 432 ± 15 MeV, which is consistent with the usual estimate (from 
extracting the string tension from the heavy quark potential) of about 440 
MeV. Using the total average they get m(0++) = 1572 ± 65 ± 55 MeV where 
the first error is statistical and the second comes from the scale. 

The IBM group, on the other hand, notices that rupa and m^a scale 
asymptotically, use the phi mass to predict quenched t^usi then extrapolate 
m(0++)/A = A + B(aA)2. They get 1740(41) MeV from their data, when they 
analyze UKQCD data, they get 1625(94) MeV, and when they combine the 
data sets, they get 1707(64) MeV. 

A neutral reporter could get hurt here. It seems to me that the lattice 
prediction for the scalar glueball is 1600 ± 100 MeV, and that there are two 
experimental candidates for it. 

Masses are not the end of the story. The IBM group has done two interest- 
ing recent calculations related to glueballs, which strengthen their claim that 
the /j(1710) is the gluebaU. 

The first one of them0 was actually responsible for the press release. It 
is a calculation of the decay width of the glueball into pairs of pseudoscalars. 
This is done by computing an unamputated three point function on the lattice, 
with an assumed form for the vertex, whose magnitude is fitted. The result 
is shown in Fig. The octagons are the results of the simulation and the 
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Figure 26: Scalar glueball decay couplings from Ref. 53. 



diamonds show interpolations in the quark mass. The "experimental" points 
(squares) are from a partial wave analysis of isoscalar scalar resonances by 
Longacre and LindenbaumE3. 

The response of a member of the other side is that the slope of the straight 
line that one would put through the three experimental points is barely, if at 
all, compatible with the slope of the theoretical points. Since they argue 
theoretically for a straight line, the comparison of such slopes is a valid one. 

If one of the experimental states is not a glueball, it is likely to be a ^Pq 
orbital excitation of quarks. Weingarten and LeeEj are computing the mass 
of this state on the lattice and argue that it is lighter than 1700 MeV; in their 
picture the /o(1500) is an ss state. I have now saicLmore than I know and will 
just refer you to recent discussions of the questionEJ. 

Both groups predict that the 2++ glueball is at about 2300 MeV. 

CarL"improved actions" help the situation? Recently, Peardon and Morn- 
ingstarc3 implemented a clever method for beating the exponential signal-to- 
noise ratio: make the lattice spacing smaller in the time direction than in the 
space direction. Then the signal, which falls like exp(— mat-Lt) after Lt lattice 
spacings, dies more slowly because at is reduced. Their picture of the glue- 
ball mass vs tq. is shown in Fig. ^ They are using the tadpole-improved 
Liischer-Weisz action. The pessimist notes the prominent dip in the middle 
of the curve; this action also has a lattice-artifact transition (somewhere); the 
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Figure 27: Glucball mass vs ro from Ref. 49, including large lattice spacing data. 

optimist notes that the dip is much smaller than for the Wilson action and 
then the pessimist notes that there is no Wilson action data at large lattice 
spacing to compare. I think the jury is still out. 

7.2 a,{Mz) 

For some time now there have been claims Ihat physics at the Z pole hints at a 
possible breakdown in the standard modelEa. A key question in the discussion 
is whether or not the value of aj^jg inferred from the decay width of the Z 
is anomalously high relative to other determinations of the strong coupling 
(when run to the Z pole, usually). 

The most recent analysis I am aware of is due to Erler and Langacker 
Ej. Currently, q,^!!^*'"^?"^ = 0.123(4) (2) (1) for the standard model Higgs mass 
range, where the first/ second/ third uncertainty is from inputs/ Higgs mass/ 
estimate of terms. The central Higgs mass is assumed to be 300 GeV, and 
the second error is for Mh = 1000 GcV (+), 60 GeV (-). For the SUSY Higgs 
mass range (60-150 GeV), one has the lower value ctj^ — .121(4)(+1 — 0)(1). 
A global fit to all data gives 0.121(4)(1). Hinchcliffe in the same compilation 
quotes a global average of 0.118(3). 

The lattice can contribute to this question by predicting ctjjg from low 
energy physics. The basic idea is simple: The lattice is a (peculiar) UV cutoff. 
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Figure 28: Survey of a-jjoiMz) from Ref. 61. 



A lattice mass /i = Ma plus an experimental mass M give a lattice spacing 
a = [ijM in fm. If one can measure some quantity related to at a scale 
Q ~ 1/a, one can then run the coupling constant out to the Z. 

ae best (recent) lattice number, from Shigemitsu's Lattice 96 summary 
. is 

a^(Z) = 0.1159(19)(13)(19) (135) 

where the first error includes both statistics and estimates of discretization 
errors, the second is due to uncertainties from the dynamical quark mass, and 
the third is from conversions of conventions. The lattice number is about one 
standard deviation below the pure Z-physics number. Lattice results are com- 
pared to other^xecent determinations of Q!p-g(Z) in Fig. |2^, a figure provided 
by P. Burrows EJ. 

Two ways of calculating_Qjs(M2) from lattice have been proposed: The first 
is the "small loop method" Ed. This method uses the "improved perturbation 
theory" described in Chapter 3: One assumes that a version of perturbation 
theory can describe the behavior of short distance objects on the lattice: in 
particular, that the plaquette can be used to define av(g = 3.41/a). With 
typical lattice spacings now in use, this gives the coupling at a momentum 
Qo = 8 — 10 GeV. One then converts the coupling-to oc-jj-g and runs out to the 
Z using the (published) three-loop beta functionE3. 

Usually, the lattice spacing is determined from the mass splittings of heavy 



51 



states. This is done because the mass differences between physical heavy 
quark states are nearly independent of the quark mass- for example, the S-P 
mass splitting of the -0 family is about 460 MeV, and it is about 440 MeV for 
the T. A second reason is that the mass splitting is believed to be much less 
sensitive to sea quark effects than light quark observables, and one can estimate 
the effects of sea quarks through simple potential models. The uncertainty in 
the lattice spacing is three to five per cent, but systematic effects are much 
greater (as we will see below). 



The coupling constant comes from Eq. 125. The plaquette can be mea- 
sured to exquisite accuracy (0.01 per cent is not atypical) and so the couphng 
constant is known essentially without error. However, the scale of the coupling 
is uncertain (due to the lattice spacing). 

The next problem is getting from lattice simulations, which arc done with 
Uf — (quenched) or — 2 (but unphysical sea quark masses) to the real 
world of nj — 3. Before simulations with dynamical fermions were available, 
the translation was done by running down in Q to a "typical gluonic scale" for 
the psi or the upsilon (a few hundred MeV) and then matching the coupling 
to the three- flavor coupling (in the spirit of effective field theories). Now we 
have simulations at Uf —2. Recall that in lowest order 



1 

as ^ 4tt ' 



{—^V^^ (136) 



One measures l/a^ in two simulations, one quenched, the other at ~ 2, runs 
one measurement in Q to the Q of the other, then extrapolates 1/a linearly in 
to = 3. Then one can convert to MS and run away. 

Pictures like Fig. ^are not very useful when one wants to get a feel for the 
errors inherent in the lattice calculation. Instead, let's run our expectations for 
OLs{Mz) down to the scale where the lattice simulations are done, and compare. 
Fig. does that. The squares are the results of simulations of charmed quarks 
and the octagons are from bottom quarks, both with n/ = 0. The crosses and 
diamond are Uf — 2 bottom and charm results. (The bursts show upsilon data 
when the 1S-2S mass difference gives a lattice spacing.) Note the horizontal 
error bars on the lattice data. Finally, the predicted = 3 coupling ap is 
shown in the fancy squares, with error bars now rotated because the convention 
is to quote an error in as. The lower three lines in the picture (from top to 
bottom) are ajjg{Mz) — 0.118, 0.123, and 0.128 run down and converted to 
the lattice prescription. 

The two top lines are predictions for how quenched a should run. 

Now for the bad news. All of the Uf — 2 data shown here were actually run 
on the same set of configurations. The bare couplings are the same, but the 
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Figure 29: Survey of a-jg^{Q) at the scale where lattice simulations are actually done. 

lattice spacings came out different. What is happening is that we are taking 
calculations at some lattice spacing and inferring a continuum numbers from 
them, but the lattice predictions have scale violations which are different. (The 
T calculations use nonrelativistic quarks, the ip calculations use heavy Wilson 
quarks.) Notice also that the bottom and charm quenched lattice spacings 
are different. This discrepancy is thought to be a failure of the quenched 
approximation: the characteristic momentum scale for binding in the tp and 
T are different, and because nf is not the real world value a runs incorrectly 
between the two scales. Said differently, in the quenched approximation, the 
spectrum of heavy quark bound states is different from the real world. 

There is a second method of determining a running coupling constant 
which actually allows one to see the running over a large range of scales. It 
goes by the name of the "Schrodinger functional," (referring to the fact that 
the authors study QCD in a little box with specified boundary conditions) but 
"coupling determined by varying the box size" would be a more descriptive 
title. It has been applied to quenched QCD but has not yet been extended to 
11 QCD. I will describe the method in the context of the d — 2 sigma model 
rather than QCDES. 
The idea is that coupling constants can be defined through the response 
of a system to boundary conditions. For example, if one dimension of a c? = 2 
0{n) sigma model is compactified with size L, the mass gap M{L) (defined 
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Table 1: Data for the running coupling constant in pure SU(2) gauge theory. 



u 


a{2,u) 


2.037 


2.45(4) 


2.380 


2.84(6) 


2.840 


3.54(8) 


3.55 


4.76 



through the transfer matrix along the other dimension) is related to the cou- 
pling through 

g^{L) = 2M{L)L/{n-l). (137) 

The coefficients insure that the answer reduces to the lowest-order result. Now 
we take take the beta function 

/3(ff)=i^, (138) 

stretch the length of the compact dimension from Lq to sLq, and compute the 
change in the coupling 

where the "step scaling function" a{s, u = g^{L)) = g'^{sL) is the new coupling 
constant. 

Now the idea is to do simulations with the same bare coupling on systems 
of size Lo and sLq and, by measuring u and a{s, u), to see how the new coupling 
depends on the original one. An example of how this works is shown in Table 
0. Each horizontal pair are the coupling and its step function (for s = 2) . The 
output coupling on the first line is used as the input coupling on the second 
line (or rather, since these couplings are derived, not bare couplings, one must 
interpolate to begin with the old output coupling as the new input coupling). 
When the matching is not perfect, perturbation theory is used for the small 
amount of running which is required. After four steps the physical length scale 
has increased by a factor of 2^ = 16. 

Now we can unfold the couplings. We measure distances in terms of the 
largest L in the simulation and show the coupling constant in Table ^j. 

We can introduce a scale in GeV by measuring something physical with the 
bare parameters corresponding to one of the L's, say, for the gauge theory, the 
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Table 2: The Schrodinger functional running coupling constant in pure SU(2) gauge theory. 



1 -^max 


9\L) 


1.0 


4.765 


0.500(23) 


3.55 


0.249(19) 


2.840 


0.124(13) 


2.380 


0.070(8) 


2.037 



Sommer variable tq (it should be obvious why the authors use this observable) 
at L„iax- An example of a running coupling constant via this prescription is 
shown in Fig. ^ This is a coupling constant in a particular prescription; 
when it is small enough it could be matched to any other prescription using 
perturbation theory. 

All the simulations are done in small lattice volumes. I have left out 
many technical details, of course. Please note that this calculation checks 
perturbation theory; it does not use it overtlvUjapart from the small amount 
of running I described). I think it is beautiful E£|. 

8 Conclusions 

There it is! I now appreciate what the "string lecturer" goes through when 
he arrives at TASI's focussed on phenomenology. I have left out a lot. For 
example, the lattice is a fertile source of numbers for comparisons with the 
standard model, where what the experimentalist measures is often a combi- 
nation of a "fundamental parameter" of the standard model times a hadronic 
matrix element. The lattice and lattice QCD in particular are going through a 
lot of changes right now, and perhaps for the better. I hope that some of the 
methods I have described might be useful as a starting point when you find 
yourselves working in the strong-coupling sector of your favorite theory. 
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